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This paper investigates whether observations of seismic anisotropy are 
compatible with a cubic structure of the inner core Fe alloy. We assume 
that anisotropy is the result of plastic deformation within a large scale flow 
induced by preferred growth at the inner core equator. Based on elastic 
moduli from the literature, bcc- or fcc-Fe produce seismic anisotropy well 
below seismic observations (< 0.4%). A Monte-Carlo approach allows us to 
generalize this result to any form of elastic anisotropy in a cubic system. 
Within our model, inner core global anisotropy is not compatible with a 
cubic structure of Fe alloy. Hence, if the inner core material is indeed cu¬ 
bic, large scale coherent anisotropic structures, incompatible with plastic 
deformation induced by large scale flow, must be present. 
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1. Introduction 

Seismic observations provide strong evidence that the Earth’s inner core is anisotropic, 
with larger velocity in the polar than in the equatorial direction [see Deuss, 2014, for 
recent review]. This anisotropy is observed both using short wavelength body waves 
differential travel times [e.g. Irving and Deuss, 2011a, and references therein], and long 
period normal modes [e.g. Irving and Deuss, 2011b, and references therein]. Over the 
years, this observation was refined with further evidences for both hemispherical [e.g. 
Tanaka, 2012, and references therein] and radial variations of the anisotropy amplitude 
and geometry, with an almost isotropic layer at the top of the inner core surrounding a 
more anisotropic region [Souriau, 2003] and possibly an innermost inner core with different 
properties [e.g. Lythgoe et al, 2014, and references therein]. 

This anisotropic structure is likely due to an alignment of anisotropic Fe-alloys acquired 
either during solidification [e.g. Bergman, 1997] or by deformation afterwards [Yoshida 
et al. [1996], and Deguen [2012] for a recent review], although other hypothesis have been 
proposed, such as the presence of melt inclusions in the solid inner core [Singh et al, 2000], 
The stable form of inner core Fe-alloy remains a matter of debate: hexagonal close-packed 
(hep), body-centered cubic(bcc) and face-centered cubic (fee) have been proposed in the 
literature. Experiments imply that the stable form for pure Fe at inner core conditions 
is hep [e.g. Tateno et al, 2010] but hrst-principles calculations also proposed various 
structures [bcc, Belonoshko et al. [2008], fee, Mikhaylushkin et al. [2007], or hep, Modak 
et al. [2007]]. However, differences in the free energy of the various Fe polymorphs are 
very close to each other [Bouchet et al, 2013], Moreover, the effect of impurity elements 
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such as Ni or Si is still not well constrained [see Morard et ai, 2014, for a recent review]. 
As such, a cubic structure for Fe at inner core conditions remains a working hypothesis. 

Recently, we developed a framework for scaling microscopic properties such as single¬ 
crystal elasticity to the macroscopic observations of anisotropy [Lincot et al, 2014]. The 
model allows us to build a synthetic inner core with a given choice of crystal structnre 
(previously hep), elastic moduli, deformation mechanisms, and geodynamical model. In 
a second stage, we simulate the observations of body wave differential travel times that 
can then be compared with seismic observations. 

Here, we investigate the effect of a cubic structure for Fe in such model and demonstrate 
that, within our framework, a cubic structure for Fe at inner conditions fails to reprodnee 
observations of seismic anisotropy in the inner core. 

2. Methods 

As in Lincot et al. [2014], our geodynamical model is based on that of Yoshida et al. 
[1996] (Fig. 1). The model assumes that geostrophic convection in the outer core induces 
faster crystallization of the inner core in its equatorial belt. The resulting topography 
is continuously relaxed by a quadrupolar flow, generating a plastic deformation (model 
Yos-Ran, Fig. 1). We also use the extension of [Deguen and Cardin, 2009] that accounts 
for a stable density stratification during inner core formation. Such model localizes the 
motions in the outer portion of the inner core, allowing mnch larger deformations (model 
Yos-Strat, Fig. 1). 

Onr choice of geodynamical model is driven by the objective of producing lattice pre¬ 
ferred orientations (LPO) aligned with the Earth’s rotation axis, leading to a large North- 
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South seismic anisotropy. To this end, quadrupolar models, such as Yos-Ran and Yos- 
Strat, are most efficient. Other models with smaller scale inner core deformation can 
result in strong LPO at the local scale and, hence, a large local anisotropy. However, such 
complex velocity structures within the inner core are averaged over the path of a seismic 
ray and will fail to produce global ’’average” anisotropy at the inner core scale [Lincot 
et aL, 2014]. 

Our study focuses on anisotropy acquired by deformation and does not account for 
anisotropy acquired during crystallization. However, crystallization textures in cubic ma¬ 
terials are complex and, hence, will probably not lead to any significant global scale 
anisotropy. Similarly, non axisymmetric models, such as those involving thermal transla¬ 
tion of the inner core, are out of the scope of this study, which focuses on models preserving 
a symmetry around the axis of rotation of the Earth. 

As in Deguen et al. [2011] and Lincot et al. [2014], we compute the position and defor¬ 
mation for 100 markers introduced at the top of the inner core boundary (ICB) during 
inner core growth (Fig. 2). Texture along these markers are calculated for a 3000 grains 
aggregate of cubic Fe using the Los Alamos viscoplastic self-consistent (VPSC) code of 
Lebensohn and Tome [1993], assuming a random crystallization texture and dominant slip 
along {111}(110) and {110}(111) for the fee and bcc structure, respectively [Kocks et al, 
1998]. Computed textures for a present day in core are presented in Fig. 2. 

A each point of the grid, the local elastic tensor of the polycristal is then calculated 
by weighting the single crystal elastic moduli with the aggregate texture under the Hill 
approximation (Fig. 4). Here, we use single-crystal elastic moduli from first-principles 
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calculations (Fig. 3), with a 3.5% and 13.3% single crystal P-wave anisotropy for the bcc 
and fee phase, respectively [Vocadlo, 2007; Vocadlo et ai, 2008]. In both cases, velocities 
are minimal and maximal along the (100) and (111) directions, respectively. Due to the 
cubic symmetry, velocity distributions in the single-crystal are complex and display 14 
extremas, unlike hep materials for which 3 extrema are found [Lincot et ai, 2014] . 

Velocity distributions, such as those presented in Fig. 4, can be difficult to compare 
with actual observations of anisotropy. Hence, following the procedure of Lincot et al. 
[2014], we generate more than 100000 synthetic seismic rays to probe the whole inner 
core and simulate seismic measurements. For each segment of the ray, we evaluate the 
slowness of the material by solving the Christoffel equation with the local elastic tensor 
deduced by linear interpolation of the elastic tensor of the aggregates located at the four 
nearest grid point of reference. For each ray, we estimate the normalized seismic travel 
times residual 

s - sO 

5t/t = — (1) 

where s is the simulated slowness of the seismic ray, and is the slowness of that same 
ray for an homogeneous and fully isotropic inner core. 

Our model allows for a detailed analysis of anisotropy, including the depth and orien¬ 
tation dependence of the travel time residuals (Fig. 5). Those, however, can be difficult 
to compare with “global” scale anisotropy reported in seismology publications. Assum¬ 
ing cylindrical symmetry of the inner core, Creager [1992] proposed to fit travel times 
residuals using 

St/t = a + bcos^ ( +Ceos‘S (2) 
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where a, b and c are adjustable parameters and ( is the angle between the ray and the 
Earth rotation axis. A measure of the global inner core anisotropy typically reported in 
the literature is the quantity b + c [e.g. Irving and Deuss, 2011a]. 

3. Results for published elastic moduli 

Without stratification (Yos-Ran, Figs. 2a,b), textures in the outer part of the inner 
core are nearly random while nearly axisymmetric textures develop in the central portion. 
This pattern is the result of the relatively slow deformation in the model of Yoshida et al. 
[1996], for which 100% deformation is reached over the age of the inner core. Pole figures 
for bcc (Fig. 2a) and fee (Fig. 2b) are inverted, with the 100 axes of the bcc phase 
preferentially aligned in the equatorial plane while, for fcc-Fe, they tend to lie in the 
North-South direction. This inversion of textures between bcc and fee structures in well 
known in materials science and is due to an exchange of slip plane and direction between 
both structures. 

With stratification (Yos-Strat, Figs. 2c,d), strong textures are mostly concentrated at 
the top of the inner core, in a superficial layer where deformation is the largest [Deguen 
and Cardin, 2009]. In this superficial layer, the orientation of the texture varies strongly 
from mid latitude, where the horizontal stress is maximum (simple shear), to polar or 
equatorial regions for which deformation is similar to that of a corner fiow. In opposition 
to model Yos-Ran, most of the LPO developed in Yos-Strat are not axisymmetric and 
heterogeneous, both in direction and in intensity in the overall inner core. 

Using the elastic moduli of Fig. 3, anisotropy in P-wave velocity is less then 0.1% in 
model Yos-Ran, whatever the location inside the inner core (Figs. 4a,b). This also applies 


DRAFT 


February 5, 2015, l:57ain 


DRAFT 



X- 8 


LINCOT ET AL.: INNER CORE MADE OF CUBIC IRON? 


to the central portion of the inner core where the LPO is the strongest. This isotropy at 
the scale of the polycrystal is the result of the multiple symmetries for velocities in the 
single-crystal (Fig. 3) that average out at the scale of the aggregate. Note that this is an 
upper estimate that neglects the solid body rotation of the inner core that could occur 
during its formation. 

The addition of stratified growth (model Yos-Strat, Figs. 4c,d) induces strong 
anisotropies in a limited region, at the base of the superficial layer at mid latitude. 
Anisotropy is more pronounced for the fee phase (Fig. 4d) reaching values of up to 
10% locally, close to that of the single crystal (13.3%, Fig. 3b). Anisotropy for an inner 
core composed of bcc-Fe with the elastic moduli of Fig 3 remains below 1.5% (Fig. 4c). 
In both cases, one should note that local anisotropy in model Yos-Strat is complex, with 
numerous extrema that do not align with a North-South axisymmetry. 

Fig. 5 presents the travel time residuals for seismic rays crossing the models of Fig. 4. 
The Yos-Ran model, without stratification, is isotropic, with a constant average ray ve¬ 
locity, whatever the direction of propagation or depth of turning point for the ray. Hence, 
the model of Yoshida et al. [1996], combined with a cubic phase of Fe and the elastic 
moduli of Fig 3 is not consistent with observations of seismic anisotropy in the inner core. 

With the addition of stratification, global anisotropy for bcc-Fe is low (0.1%, Fig. 5c) 
and opposite to that observed in seismic studies, i.e., North-South ray paths are slower 
than those in the equatorial plane. For fcc-Fe, a larger global anisotropy of 0.4% (Fig. 5d) 
is generated, as expected from Fig. 4d. Nevertheless, note that the global anisotropy is a 
rather moderate in relation to that of the single crystal (Fig. 3) and the seismic reports of 
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over 3% global anisotropy in the literature [Irving and Deuss, 2011a]. The local maxima of 
P-wave anisotropy of over 10% observed in Fig. 4d are averaged out over the integration 
along the ray path. 

Residuals in Fig. 5d are minimal for rays with ( = 60°, in opposition to actual seismic 
observations for which residual remain fairly constants for ( > 45° [Irving and Deuss ^ 
2011a]. The (F-shape of the calculated residuals in model Yos-Strat is mainly due to the 
contribution of rays with a turning point between 244 and 488 km below the ICB (blue 
dots),that is for rays that probe just below the first stratified layer with strong anisotropy 
in Figs. 4c,d. This (F-shape also implies that North-South and equatorial rays travel with 
a relatively fast speed. For deeper rays (red and black dots in Fig. 5), the variations 
of the travel time anomalies are smaller with C, (less than 0.5%) but North-South rays 
travels faster than others. Also note that the fit of Eq. 2 is a poor representation of the 
travel time residuals and that the scattering of the data points around the fit is rather 
large (more than 1% of 5t/t, with a global anisotropy b -\- c — 0.4%). 

In summary, we find that, within our synthetic grown inner core, no cubic phase of 
iron (bcc nor fee) with the elastic moduli from first-principles calculations [Vocadlo^ 2007; 
Vocadlo et al, 2008] can generate a global anisotropy that compares with observations in 
seismology [e.g. Irving and Deuss, 2011a]. 

4. Generalisation to other elastic models 

Our conclusions from the previous section are strongly constrained by the choice of 
elastic moduli in Fig 3. This section intends to generalise our result to any cubic phase, 
whatever the choice of elastic moduli. To that end, we introduce a Monte-Carlo approach 
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to determine the effect of single-crystal elasticity on the global anisotropy at inner core 
scale. 

We select 3000 random sets of elastic moduli constrained by conditions for mechanical 
stability [Wallace, 1972] 

Cii — Ci2 > 0, Cii 2 Ci2 > 0, (7*44 > 0. (3) 

and an average bulk and shear moduli in the Hill average that match those of PREM [K 
= 1384 GPa, G = 166 GPa Dziewonski and Anderson, 1981] with 15%. These 3000 sets 
of elastic moduli sample all possibilities for single-crystal elasticity for cubic Fe at inner 
core conditions. 

Single-crystal anisotropy can be quantified using the dimensionless parameter K, as 
defined in Karki et al. [1997]: 

K = 1 

Gil 

For an isotropic crystal, K = 0. K becomes positive (negative) when P-wave velocities 
are minimal (maximal) along (100). With the elastic moduli of Fig. 3, K is positive. 

For each of the 3000 sets of elastic moduli, we repeat the procedure above to characterize 
the global inner core anisotropy for each of the 4 models, i.e. geodynamical models Yos- 
Ran and Yos-Strat, with an inner core consisting of either a bcc- or an fcc-structured 
material. Global anisotropy (6 -|- c) plotted as a function of the dimensionless parameter 
K follows a linear trend (Fig. 6): global anisotropy in the inner core increases with single 
crystal anisotropy. 

Without stratification, no cubic material can produce sufficient global anisotropy to 
match seismic observations. This result is general and does not depend on the crystal 
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structure, nor the set of elastic moduli. Stratification enhances global inner core anisotropy 
(model Yos-Strat, Fig. 6c,d) but global anisotropy does not exceed 0.5%, well below 
seismic observations. Anisotropies are larger for the case of fee structure. 

Seismic observations of inner core anisotropy are heterogeneous, depth, and geographi¬ 
cally dependent [Irving and Deuss, 2011a]. Hence, a global fit based on Eq. 2 on 100000 
random rays may not be relevant for comparing with seismic observations. Hence, for each 
of the 3000 random sets of elastic moduli. Fig. 6 also presents the standard deviation of 
all individual residuals. Using Yos-Ran model, none of the residuals will ever exceed 0.2% 
whatever the elastic moduli, for all ray paths. For the Yos-Strat model, standard devia¬ 
tions of all residnals can reach values of up to 1% for extremely anisotropic single crystal 
elastic moduli. Hence, seismic observations of above 3% anisotropy are incompatible with 
our model, whatever the ray path, whatever the elastic moduli, and whatever the choice 
of crystal structure, bcc or fee. If the inner core is cubic Fe, its seismic anisotropy cannot 
be due to plastic deformation alone. 

5. Conclusions 

We used the framework of Lincot et al. [2014] to study whether inner core anisotropy 
arises from plastic deformation of a cubic phase. Our analysis relies on a geodynamical 
model derived from that of Yoshida et al. [1996], Fe alloy either in the bcc or fee structure 
deforming along classical slip systems, and a Monte-Carlo approach for probing all possible 
sets of elastic modnli. 

Within onr assnmptions, no cubic phase of Fe can produce a global inner core anisotropy 
larger than 0.5%, far below seismic observations. This observation results from the mul- 


DRAFT 


February 5, 2015, l:57ain 


DRAFT 



X - 12 


LINCOT ET AL.: INNER CORE MADE OF CUBIC IRON? 


tiple symmetries over multiple scales involved: symmetries in single-crystal P-wave veloc¬ 
ities and multiple equivalent slip systems combined with large inner core portions probed 
with each seismic ray. 

Hence, the observed inner core anisotropy does not arise from deformation-induced plas¬ 
tic deformation of a cnbic phase of an Fe alloy. If the inner core material is indeed cubic, 
other mechanisms for generating anisotropy are required and preserved, such as solidih- 
cation texturing, grain growths or melt inclusions. In any case, anisotropic structures in 
the inner core should be coherent over large scales for a cubic material to be compatible 
with observations of seismic anisotropy. 
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Figure 1. Inner core geodynamical models used in this study. In Yos-Ran, a quadripolar 
flow results from the preferred equatorial growth of the inner core [Yoshida et ai, 1996]. 
Yos-Strat adds the contribution of a stable density stratification, resulting in layered flow 
pattern [Deguen and Cardin, 2009]. 
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r/r,c r/r,c 


Figure 2. 100 pole figures representing present-day textures in Fe aggregates at various 

locations in the inner core for models Yos-Ran {a,h) and Yos-Strat {c,d) and either bcc-Fe 
{a,c) or fcc-Fe {b,d). The vertical z axis is the geographical South-North axis while x lies 
in the equatorial plane. Both geodynamical models are symmetric about x and axisym- 
metric around z. Blue lines are the trajectories of the polycrystalline aggregates after 
crystallization at the ICB. Contours for pole figures in multiples of a random distribution 
(m.r.d.). 
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bcc-Fe fcc-Fe 

[Vocadio et a!., 2007] [Vocadio et al., 2008] 



<010> Max anisotropy: 3.5% Max anisotropy: 13.3% 



Figure 3. P-wave velocity (in km/s) in bcc-Fe (left) and fcc-Fe (right) single crystal 
at inner core conditions (p = 13155 kg/m^, T = 5500 K) corresponding to single crystal 
elastic moduli (in GPa) of Cn = 1505, C 12 = 1160 and C 44 = 256 for bcc-Fe [Vocadio, 
2007] and Cn = 1397, C 12 = 1247 and C 44 = 423 for fcc-Fe [Vocadio et al, 2008]. Inlet on 
the right shows elastically equivalent directions in a cubic structure such as ( 100 ), ( 010 ), 


( 001 ), ( 100 ) and ( 010 ). 


Figure 4. P-wave velocity (in km/s) in the present day inner core for geodynamical 
models Yos-Ran {a,b) and Yos-Stmt {c,d) and either bcc-Fe (a,c) or fcc-Fe {b,d). Cal¬ 
culations are based on the textures of Fig. 2 and the single crystal elastic moduli of 
Fig 3. 
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Figure 5. Travel time residuals (dots) vs. angle between the ray and the Earth’s 
rotation axis. Colors indicate the depth of turning point of the ray. Thick purple line is 


the fit of Eq. 2 to the data. Numbers in the insets indicate the global anisotropy for each 
model. Results plotted for geodynamical models Yos-Ran {a,h) and Yos-Stmt {a,b) and 
either bcc-Ee {a,c) or fcc-Ee {b,d). Calculations based on the velocities of Eig. 4 and the 
single crystal elastic moduli of Fig 3. 
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Figure 6. Global inner core anisotropy vs. single crystal anisotropy for geodynamical 
models Yos-Ran {a,b) and Yos-Stmt {c,d) and either bcc-Fe {a,c) or fcc-Fe {b,d). For 
each figure, calculations were performed on 3000 random sets of elastic moduli covering 
all possible values for the inner core (solid black dots, see text). The error range (red area) 
is an estimate of the spread of residuals for individual rays. White-filled diamond and 
yellow-filled squares are the results of calculations based on the elastic moduli of Fig 3. 
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